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Abstract. Supersymmetric quantum mechanics is a powerful tool for generating exactly 
solvable potentials departing from a given initial one. If applied to the harmonic oscillator, 
a family of Hamiltonians ruled by polynomial Heisenberg algebras is obtained. In this paper 
it will be shown that the SUSY partner Hamiltonians of the harmonic oscillator can produce 
evolution loops. The corresponding geometric phases will be as well studied. 
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1 Introduction 

In the last decades there has been a growing interest in studying evolution loops (EL), which 
are circular dynamical processes such that the evolution operator of the system becomes the 
identity at a certain time [30, 33, 37, 39, 40, 44, 45, 48, 50, 51, 53, 54]. They represent a natural 
generalization to what happens for the harmonic oscillator. Their importance rests on the fact 
that the EL are quite sensitive to external perturbations, so they are a good starting point to 
implement the dynamical manipulation for approximating an arbitrary unitary operator [48] (see 
also [29]). 

On the other hand, the polynomial Heisenberg algebras (PHA) are deformations of the 
Heisenberg-Weyl algebra in which the commutators of the Hamiltonian H with the annihi- 
lation L~ and creation L + operators are standard but the commutator between the last two 
turns out to be a polynomial in H [1, 3, 6, 9, 13, 14, 17, 18, 27, 35, 38, 46, 47, 49, 60]. In order 
to characterize the spectrum of H, Sp(H), one needs to determine the eigenstates of H which 
are extremal (annihilated by L~) and have physical interpretation: thus, Sp(-ff) is composed of 
several independent ladders (either of finite or infinite lengths) departing from those extremal 
states (for alternative deformations of the Heisenberg- Weyl algebra see e.g. [23, 25, 42, 57]). 

In addition, nowadays it is widely accepted that supersymmetric quantum mechanics (SUSY 
QM) is the simplest technique for generating new Hamiltonians H^ departing from a given 
initial one Hq (for recent books and review articles see [5, 10, 12, 21, 26, 31, 34, 41, 55, 59]). 
After applying the method, it turns out that Sp(-fffc) differs little from Sp(-ffn) (the differences 
rely in a finite number of levels). Moreover, the algebraic structure of H^ can be obtained 
straightforwardly of the corresponding algebra of Hq. In particular, the SUSY partners of the 
harmonic oscillator Hamiltonian turn out to be ruled by polynomial Heisenberg algebras [35, 49]: 
this is the simplest way for realizing in a non-trivial way such non-linear algebras. 
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Figure 1. A system which performs an evolution loop at t = r when perturbed, in general, will deviate 
from this cyclic dynamical process. 

In this paper I would like to explore in detail the possibility that the SUSY partners Hamil- 
tonians of the harmonic oscillator can have EL, as it happens for the initial system. If the 
answer becomes positive, it is natural to evaluate then the geometric phase, associated to an 
arbitrary initial state which is cyclic [2, 4, 15, 16, 19, 30, 39, 40, 44, 45, 56, 58]. It is worth 
to notice that some partial results along this way have been derived previously [30]. However, 
they were found just for the Abraham-Moses family of potentials isospectral to the harmonic 
oscillator [49] . Here, we are going to generalize these results for a SUSY transformation of order 
k > 1 so that H$ and are not necessarily isospectral [35]. 

The paper has been organized as follows. In the next section we shall introduce briefly 
the evolution loops, with a discussion about the geometric phases which naturally will arise 
for systems satisfying such operator identity. In Section 3 the polynomial Heisenberg algebras 
will be presented in general, while in Section 4 the same shall be done for supersymmetric 
quantum mechanics. The SUSY partners of the harmonic oscillator will be derived at Section 5, 
including the analysis of their connections with polynomial Heisenberg algebras, evolution loops 
and associated geometric phases. Our conclusions shall be presented at Section 6. 



2 Evolution loops 

At operator level, the dynamics of a quantum system is determined by its evolution opera- 
tor U(t), which satisfies: 

= -iH(t)U{t), U(0) = I, 



dt 



was 



where H(t) is the system Hamiltonian, / is the identity operator, U(t) is unitary. As it 

pointed out previously, we are interested in studying system having evolution loops, i.e., circular 
dynamical processes such that U (t) becomes the identity (up to a phase factor) at a certain time, 

U(r) = e*I, (1) 

with r > being the loop period, tp E M [30, 48, 51]. The simplest system with an evolution 
loop of period r = 2ir is the harmonic oscillator since 



oo 



n=0 n=0 

oo 

E/(T) = e-^X>n)<V>n|=-/, 



n=0 
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Figure 2. Representation of the cyclic state \ip(t)) £ H induced by Schrodinger equation, its closed 
shadow Tr(t) £ V and the corresponding horizontal lifting \ip(t)) £ H. The holonomy of this lifting is the 
geometric phase factor . 

where we have employed that H\ip n ) = E n \ip n ) with E n = n + 1/2, n = 0, 1, . . . (we are using 
natural units such that H = m = co = 1). The EL are important since they can be used as 
a starting point to implement the dynamical manipulation in order to approximate any unitary 
operator [48, 50, 51]. In fact, there is a prescription for implementing this kind of manipulation, 
which was introduced some years ago [48]: first of all the system has to be placed in an EL, 
i.e., Uq(t) = /; by perturbing then the EL, the small deviations of this dynamical process such 
that U(t) / will eventually approximate an arbitrary unitary operator (see an illustration in 



It is worth to notice that the EL have been mainly studied for systems ruled by time- 
dependent Hamiltonians, either in one or several dimensions [37, 45, 48, 50, 51, 53, 54] or 
for purely spin systems [39, 40, 44]. However, there are several works where the evolution loops 
are produced by time-independent Hamiltonians [30, 33]. In particular, an interesting physical 
system of such a type consists of a charged particle inside an ideal Penning trap [33] . 

Since it is an operator relationship (compare with [28]), the requirement of equation (1) is 
quite strong: it implies that for a system having an EL any £ T~L, taken as an initial condition, 
becomes cyclic with period r: 



Thus, it is natural to ask if the global phase (p has a geometric component. As it was noticed 
by Aharonov and Anandan [2], it turns out that a geometric phase f3 can be associated to any 
cyclic evolution \ip(t)) satisfying equation (2): 



J o 

In particular, if the system is ruled by a time-independent Hamiltonian H(t) = H the previous 
expression becomes simpler [30]: 



Let us note that j3 describes a global "curvature effect" arising on the space of physical states 
of the system, which is the projective space V formed by the rays or the density operators \ip){ip\ 
instead of the Hilbert space T~L [2, 16, 30, 39, 40]. Due to this curvature, the horizontal lifting 
(parallel transport) of the closed trajectory ir(t) = |^(t))(V , (*)| £ V leads to a trajectory \ip(t)} 
which is, in general, open on T~L. The holonomy of this lifting is the Aharonov-Anandan geometric 
phase factor e l/3 (see Fig. 2). 



Fig. 1). 



\^r)) = U(r)m = e^). 




/3 = ^ + t(V>|H|V>). 



(3) 



4 



D.J. Fernandez 



: • : • Sj+i-i 

*l— — — = 



(a) 



(b) 



Figure 3. Spectrum of the Hamiltonians ruled by polynomial Heisenberg algebras: (a) it is composed 
of s infinite ladders starting from £j, j = 1, . . . , s; (b) the j-th ladder becomes finite since equation (7) 
is satisfied. 



3 Polynomial Heisenberg algebras 

The polynomial Heisenberg algebras are deformations of the Heisenberg- Weyl algebra of kind 
[17, 35]: 

[H,L ± ] = ±L ± , (4) 
[L~,L+] = Q m +i(H + 1) - Q m +i(H) = P m (H), (5) 

where 

m+l 

Q m+ i(H) = L+L- = 11 (H - Si) (6) 

i=i 

is a (m + l)-th order polynomial in H which implies that P m {H) is a polynomial of order m-th 
in H. A simple way of realizing the algebra of equations (4)~(6) is to suppose that H has the 
standard Schrodinger form, 

while are (to + l)-th order differential operators. 

Note that Sp(-ff) depends on the number of eigenstates of H belonging to the kernel of L~ 
which have physical meaning. If s of them are physically acceptable and satisfy 



L ifj £i = 0, Hi/j £i = , i = l,...,s, 

thus, Sp(H) turns out to be composed of s independent infinite ladders, each one of them 
starting from tpg i , i = 1, . . . , s (see Fig. 3a). 

On the other hand, for the j-th ladder which starts from i/j£ it could happen that 

(L+y- 1 ^ + 0, (L+)^ £ . = 0, (7) 

for some integer I. In this case it turns out that [35] 

£ n = £j+l, n € {s + l,...,k}, 

which means that the j-th ladder starts from the eigenvalue £j and ends at £j + I — 1, i.e., it is 
a finite ladder of length I, with I steps (see Fig. 3b). 
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4 Supersymmetric quantum mechanics 

Let us consider now the following chain of intertwining relationships: 

HiAf = AfHi-i, (8) 

i = l,...,k, (9) 



1 ~ V2 



d . . 
T-r +ai{x,€i) 



where 



1 d 



2 



Hi = ~2d^ + Vi{X) ' i = >---: k - (10) 

By plugging the expressions (9), (10) into equation (8), it turns out that the following must be 
satisfied: 

a'^x, £i) + af(x, €i) = 2[Vi-i(x) - (11) 
V i (x) = V i - 1 (x)-a' i (x,e i ). (12) 

Suppose now that Vi-i(x) is known; then Vi{x) becomes determined (see equation (12)) if the 
solution a.i(x, ei) of the i-th Riccati equation (11) associated to ei can be found. The key point 
in this treatment is to realize that there is a simple finite difference formula allowing to find 
algebraically cti(x, ej) in terms of two solutions of the (i — l)-th Riccati equation, associated to 
the factorization energies ej_i, e% [36, 52]: 

/ \ / \ 2(ej_i - €i) 

oti{x,ei) = -ai-i(x,ei-i) 



oti-x{x,ei-x) - Oi-\{x,ei)' 

By iterating down this equation it turns out that, at the end, cti(x, €j) can be expressed in terms 
of the i solutions 

ai(x,ej) + af(x,ej) = 2[V (x) - ej], j = 1, . . . ,i, 

of the initial Riccati equation, or in terms of the corresponding solutions of the Schrodinger 
equation, 

--u" + V (x)uj = ejUj, j = l,...,i, (13) 

where ot\ (x, ej ) = u'J Uj . 

In order to connect the previous technique and supersymmetric quantum mechanics, let us 
realize now the standard SUSY algebra with two generators 

[Q,,H ss ]=0, {Qi,Qy} = %H ss , ij = l,2, 

in the following way [8, 7, 11, 20, 22, 31, 32, 34, 43, 55] 

where 

B+B k = {H k -e 1 )---{H k -e k ), B k B+ = (H - e x ) ■ ■ ■ (H - e k ), (14) 
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Hq and H k being the initial and final Hamiltonians, intertwined by fc-th order differential 
intertwining operators, namely, 

H k B+ = B+H , H B k = B k H k , B+ = A+ ■ ■ ■ A+ , B k = A[ ■ ■ ■ . (15) 

The initial and final potentials Vq, V k , are interrelated by: 

k 

V k {x) = V (x) - e i) = Vo(x) - {HW(u l7 u k )}}", 

i=i 

where W(u±, . . . , u k ) is the Wronskian of the k Schrodinger seed solutions u±, . . . , u k . Let us note 
a certain resemblance of the fc-th order SUSY QM presented here with the method of fractional 
supersymmetry discussed elsewhere [24]. 

The previous technique has been employed successfully to generate new solvable poten- 
tials V k (x) departing from a given initial one Vq{x) for several interesting physical systems 
[10, 21, 31, 55]. Of our particular interest is the case of the harmonic oscillator [13, 17, 34, 35, 
36, 38], which is worth of an explicit discussion. 



5 Harmonic oscillator SUSY partners 

In order to implement the SUSY technique, we need to find first the general solution of the 
stationary Schrodinger equation (13) for Vo(x) = x 2 /2 and an arbitrary e, which turns out to 
be: 



u 



e 2 



1 - 2e 1 



'2'' 



; x + 2xv 



£(g) 
r(^) 



3-2e 3 



2 ,X 



(16) 



where iFi(a,b;y) is the (Kummer) confluent hypergeometric function. Let us perform now 
a non-singular k-th order SUSY transformation which creates precisely k new levels, by simplicity 
placed below the ground state energy Eq = 1/2 of the oscillator [35]. If the factorization energies 
are ordered as e k < e^-i < • • • < e± < 1/2, then the non-singular SUSY transformations arise 
for \vi\ < 1, | vi | > 1, 1 1/3 1 < 1, . . . The spectrum of the corresponding Hamiltonian reads: 



Sp(H k ) = {e k , . . . , ei, E n = n + 1/2, n = 0, 1, . . . }. 



(17) 



An illustration of a second-order SUSY partner potential of the oscillator for (ei, 62) = (—1, — §) 
and (fi, v^) = (0, 2) is shown in Fig. 4. 

Let us note that, for the Hamiltonian H k , there exists a natural pair of ladder operators: 



L7 = B+aB k , L+ = B+a + B k , 



(18) 



which are differential operators of order {2k + l)-th satisfying [17, 35, 49]: 



[H k ,L±] = ±L ± 



k ■ 



(19) 



By making use of the intertwining relationships (15) and the factorizations of equations (14), it 
is straightforward to show that: 



Q2k+i{H k ) — L^L k 



Hk 



7 i=i 



€i - 1) (H k - Ei) 
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Figure 5. Diagram representing the spectrum and the action of the ladder operators for the SUSY 
partner Hk of the harmonic oscillator Hamiltonian Hq. They are built up from the intertwining opera- 
tors Bk, Bt and the harmonic oscillator annihilation and creation operators a, a + (see equation (18)). 



This implies that the operators {Hk, L k , } generate a polynomial Heisenberg algebra of order 
(2fc)-th, which is characterized by equation (19) and the deformed commutator: 

From the analysis of the roots involved in Q2k+i{Hk), it turns out that the Sp(Hk) given in 
equation (17) can be seen as containing k + 1 ladders: k of them are one-step ladders, starting 
and ending at ej, j = l,...,k; in addition, there is an infinite one starting from | [17, 35]. 
A representation of Sp(-EZfc) and the actions of the ladder operators are given in Fig. 5. 

We have already all the elements for answering the main question we would like to pose in 
this paper: since the harmonic oscillator Hamiltonian has an evolution loop, it is natural to ask 
if its SUSY partners show as well such a kind of closed dynamical processes. First of all let us 
write down the evolution operator U(t) = e~ lHkt associated to H^: 

k oo 
j=l n=0 
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If the k factorization energies ej, j = 1, . . . , k are arbitrary, it turns out that a partial evolution 
loop is produced for r = 2ir since: 

k oo 

17(2*) = *£e-^\^.)(^.\ + e~™ £ |^)(^|. 

j=l n=0 

This means that any state belonging to the subspace generated by {\ip^),n = 0, 1, . . . }, 

oo oo 

|V(o)> = J>h/£), $> n | 2 = i, 

n=0 n=0 

becomes cyclic with period r = 2ir: 
Mr)) = e-** |V(0)). 

A straightforward calculation leads now to the associated geometric phase (see equation (3)): 

oo 

/3 = -vr + 27r(V(0)|F fe |^(0)) =2vr^n|c n | 2 . (20) 

n=l 

In particular, if c n = 5 n ^ m it turns out that 
(3 = [mod(2vr)]. 

It is interesting as well to evaluate the geometric phases associated to the coherent states which 
are eigenstates of the annihilation operator of equation (18) [35], namely, 

L^\z) = z\z), z G C. 

By expressing \z) in terms of the eigenstates of H^, 

k oo 
j=l n=0 

it turns out that 



bj — 0, c n — , c w — i. 

i/n Y[(n-€i- i)(n-ei + i) 



By iterating down the recurrence relationship for c n , it turns out that it becomes expressed in 
terms of cq. The last coefficient is fixed from the normalization condition, leading to: 



z n \^ 



n=0 



'n!nr(n + i- ei )r(n+|- ei ) 
i=i 



N(r) 



nr(i- ei )r(|- ei ) 

i=i 



\ oF2k{\ — £ i> • • • i \ — e k, | — ei, • • • , | — r 2 ) ' 
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Figure 6. Plot of the geometric phase f3 of equation (22) as a function of r (black curve). The geometric 
phase associated to the standard coherent states j3 = 2ttt 2 is as well plotted (gray curve). 



where r = \z\. By using equation (20), the associated geometric phase becomes now: 



P 



2vrH 



> 2 ^fe> 2 ^1) • • • J 2 r ) 



n(i-«xi-*) 



0-^2fc(2 ~~ £ ii • 



■ • J 2 



' 2 



■ 



■ ■ • , 2 



(21) 



Let us note that, for k = 1 and ei = —1/2, the expression of equation (21) reduces (mod (2tt)) 
to the expression of equation (27) of [30] since 



Pgcs - 2-7T = 2-7T 



0^2(1, 
0^2(1, 2; r 2 ) 



irr 



2 F 2 (2,3;r 
0^2(1, 2; r 2 )" 



On the other hand, for (ei, e 2 ) = (—1,-1) it turns out that 



1607rr 2 0^4(1,1, j,j;r 2 ) 

1^77 P/3 5 17 27. ^2^ 
1611 0-^4(2) 2 ' 10' 10 > / 



(22) 



A plot of this geometric phase as a function of r is shown in Fig. 6 (black curve) . The geometric 
phase acquired by the standard coherent states, which turns out to be j3 = 2irr 2 , is as well drawn 
(gray curve). 

Coming back to our general subject, if the factorization energies are such that 



2 rrij 



V j = 1, . . . , k, 



where L, rrij are coprime, then a global evolution loop of period r = 2M7T is obtained, namely, 



J =1 

k 



n=0 



n=0 



with M being the least common multiple of {rrij, j = 1, . . . , k}. In this case any arbitrary initial 
state 



i^(o)) = e^i<) + E c «i^)' 

j=l n=0 



j=l n=0 
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is cyclic with period r = 2Mtt: 

mr)) = e- iM -\i;(0)}- 
The associated geometric phase becomes finally: 

k oo / -. \ \ / oo k 



(3 = -Mvr + 2M7T £ e#,f + £ in + - ) |c n | 2 = 2Mtt J] n|c„| 2 - £ -^f 
\i=i n=0 / \n=l i=l J 



6 Conclusions 

In this paper it has been shown that the SUSY partners of the harmonic oscillator Hamiltonian 
realize straightforwardly the polynomial Heisenberg algebras of order 2k. As a consequence, if 
the SUSY transformation creates k new levels for H^, the corresponding spectrum can be seen 
as composed of k + 1 independent ladders, k of them being one-step ladders and an infinite one 
starting from Eq = 1/2. It has been proven also that the corresponding Hamiltonians present 
in general the so-called partial evolution loops, which induce cyclic evolutions in the subspace 
generated by the eigenstates associated to E n and, consequently, have associated geometric 
phases. This applies, in particular, to the coherent states which are eigenstates of the natural 
annihilation operator of the system, for which a general formula for the associated geometric 
phase has been derived in this paper. In particular, from this general result we have recovered 
the expression for the geometric phase which was found previously [30] for the Abraham-Moses 
family of potentials isospectral to the harmonic oscillator. Finally, we have shown that it is 
possible to produce global evolution loops by imposing restrictions on the involved factorization 
energies. The associated geometric phases induced by the last operator identity have been as 
well evaluated. 
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